Low-energy structure of the intertwining double-chain ferrimagnets 
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Motivated by the homometallic intertwining double-chain ferrimagnets A3Cus(P04)4 (A = 
Ca, Sr, Pb), we investigate the low-energy structure of their model Hamiltonian Ti — ^2 [Ji(S„-.i + 
Sn-.3) + J2(S n +i:i + Sn-i-.s)] ■ S n :2, where S „ : i stands for the Cu 2+ ion spin labeled I in the nth 
trimer unit, with particular emphasis on the range of bond alternation < J2/J1 < 1. Although 
the spin-wave theory, whether up to 0(S 1 ) or up to O(S ), claims that there exists a flat band in 
the excitation spectrum regardless of bond alternation, a perturbational treatment as well as the 
exact diagonalization of the Hamiltonian reveals its weak but nonvanishing momentum dispersion 
unless J2 = J\ or J2 = 0. Quantum Monte Carlo calculations of the static structure factor further 
convince us of the low-lying excitation mechanism, elucidating similarities and differences between 
the present system and alternating-spin linear-chain ferrimagnets. 

PACS numbers: 75.10.Jm, 75.50.Gg, 75.40.Cx 



I. INTRODUCTION 

It is a long-standing and still challenging theme in 
materials science to design molecular systems ordering 
ferromagnetically. 1 The naivest idea of ferromagnetically 
coupling nearest-neighbor magnetic centers leads to the 
highest spin multiplicity but critically depends on some 
structural parameters which are hard to handle chemi- 
cally. An alternative solution to highly magnetic ground 
states consists of aligning molecular bricks so as to ob- 
tain a nonzero resultant spin in the ground state and then 
coupling the chains again in a ferromagnetic fashion. A 
variety of quasi-one-dimensional ferrimagnets were thus 
synthesized and not a few of them have been attracting 
theoretical as well as experimental interest. 

Bimetallic chain compounds are early examples and 
among others is MnCu(pbaOH)(H 9 0) 3 (pbaOH = 
2-hydroxy-l, 3-propylenebis(oxamato) = C7H6N2O7), 2 
which retains the long-range ferromagnetic order on 
the scale of the crystal lattice. Replacing the Mn 2+ 
ions by Fe 2+ , Co 2+ , and Ni 2+ ions, Kahn and co- 
workers further synthesized a series of isomorphous 
compounds, 3 which stimulated extensive chemical ex- 
plorations of heterometallic chain magnets^ 5 - and sys- 
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FIG. 1: Cu 2+ trimeric chains in v4 3 Cu 3 (P04)4. The strongly 
coupled Cu + trimer consists of a central square planar 
Cu 2+ (1) ion (black circle) and two pyramidal Cu 2+ (2) ions 
(gray circles) bridged by oxygen ions (open circles). 
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tematic theoretical investigations of alternating-spin 
chains i^^ ' 10 ! 11 ) 12 ! 13 ! 14 In an attempt to obtain sub- 
stantially larger couplings between neighboring mag- 
netic centers and possibly attain transitions to three- 
dimensional order at higher temperatures, Caneschi et 
alJ^ made a distinct attempt to bring into interac- 
tion metal ions and stable organic radicals. The repre- 
sentative materials of general formula Mn(hfac) 2 NIT-_R 
(hfac = hcxafluoroacetylacetonate = C5H2O2F6; 
NIT-i? = nitronyl nitroxide radical = C 7 Hi 2 N 2 02-i? 
with R = CH3,C2H5,C3H 5 ,C6H 5 ) indeed exhibit anti- 
ferromagnetic intrachain interactions ranging from 200 to 
330 cm -1 . The metal-radical hybrid strategy, combined 
with fabrication of novel polyradicalsr^- yielded various 
polymerized heterospin chain compounds ! 17 : 18 

Homometallic ferrimagnetism is also realizab le) 19 ' 20 but 
its mechanism is often more subtle, essentially depend- 
ing on the structural features of the system. Coron- 
ado et al£l-£2- pioneeringly synthesized chain-structured 
compounds of such kind, M 2 (EDTA)(H 2 0) 4 • 2H 2 
(M = Ni, Co; EDTA = ethylenediamminetetraacetate = 
C10N2O8), whose ferrimagnetic behavior originates from 
the alternating g factors and is therefore faint. Ho- 
mometallic chain compounds of more pronouncedly fer- 
rimagnetic aspec t 23 : 24 ' 25 : 26 were not obtained until an- 
other decade had passed, where particular topologies 
were elaborately imposed on the intrachain exchange 
interactions. A series of compounds, M(i?-py) 2 (N 3 )2 
(M = Cu, Mn; R-py = pyridinic ligand = C 5 H 4 N-i? 
with R = CI, CH3, • • •), consists of bond-polymerized ho- 
mometallic chains, where the neighboring metal ion spins 
are bridged by versatile azido ligands and are coupled to 
each other ferromagnetically or antiferromagnetically. 

The homometallic intertwining double-chain com- 
pounds v4 3 Cu3(P0 4 )4 {A = Ca, Sr, Pb)^^^ which are 
illustrated in Fig. [TJ are topological ferrimagnets^ - in the 
strict sense. Their hybrid analogs Caa-^Sr^CuapO^ 
(0 < X < 3)2£ were also fabricated in an attempt to 
tune the antiferromagnetic bridges between the Cu(l) 



and Cu(2) sites, labeled Ji and J2, and possibly ex- 
plore how paramagnetic spins grow into bulk ferrimag- 
nets. The magnetic centers without single ion anisotropy 
and the simple crystalline structure without any organic 
ligand will contribute toward revealing intrinsic features 
of one-dimensional ferrimagnetic phenomena. Thus mo- 
tivated, various experiments have been performed on 
these copper phosphates in recent years, including high- 
field magnetization,?! specific-heat^ inelastic neutron- 
scattering j22, nuclear spin-lattice relaxation-time^ and 
electron-spin-resonance 3 ^ measurements. 

It is therefore unfortunate that theoretical investiga- 
tions of this system still stay in their early stagei 30 ' 36 i 37 
Indeed there exists a field-theoretical study 38 deserving 
special mention, but the authors restricted their argu- 
ment to the particular case of J\ = J2 taking a main 
interest in realizing organic ferromagnetism. A recent 
numerical diagonalization study 3 ^ is also a fine guide 
to this system, but the authors still devoted themselves 
to clarifying the electronic correlation effect on unsatu- 
rated ferromagnetism rather than geometrically modify- 
ing this unique bipartite lattice, starting from a model 
of the Hubbard type. An introduction of bond alterna- 
tion 5 = J2/J1 ^ 1 to this system will not only con- 
tribute toward understanding the magnetic properties of 
A 3 Cu 3 (POi)^^^^ but also illuminate the character- 
istic of the uniform point 5 = 1. We are thus led to report 
the whole excitation mechanism of homogeneous-spin 
intertwining double-chain ferrimagnets, employing both 
analytical and numerical tools. According to the spin- 
wave theory, there exist three modes of elementary exci- 
tation, two of which exhibit parallel dispersion relations, 
while the rest of which is of no dispersion, regardless 
of bond alternation. However, the exact-diagonalization 
and perturbational calculations disprove the spin-wave 
scenario that the low-lying excitation spectrum remains 
qualitatively unchanged with varying 5. Indeed there 
exist local excitations which are rigorously immobile at 
5 = 1, but they can be itinerant with 5 moving away from 
unity. Except for the two particular points 5=1 and 
(5 = 0, corresponding to a plaquette chain and decoupled 
trimers, respectively, there is no flat band in the exci- 
tation spectrum of the homogeneous- spin trimeric chain. 
We further inquire into thermal excitations based on such 
an energy spectrum. Calculating the static structure fac- 
tor as a function of temperature for an alternating-spin 
linear-chain ferrimagnet as well as for the present system, 
we show what are the universal ferrimagnetic features 
and how they vary with decreasing 5. 



II. PLAQUETTE CHAINS 

The Hamiltonian of our interest is represented as 

N 
H = Hi + H2 = y jJl(S n g + S n - 3 ) ■ S n: 2 



where S n: i symbolizes the Cu 2+ ion spin (S = |) labeled 
I in the nth trimer unit (see Fig. [T]) and the intratrimcr 
( Ji) and intertrimer (J 2 ) exchange interactions, denoted 
by Hi and Hi, respectively, are defined as < J2 < J\. 
First we take a look at the particular point of J2/J1 = 
5=1, where the model reads a plaquette chain, bearing 
some analogy with a linear chain of alternating spins 1 
and i 

Introducing bosonic operators through the Holstein- 
Primakoff transformation 
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defining their Fourier transforms as 
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a>n:U 
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with the lattice constant set equal to unity, and further 
processing them via the Bogoliubov transformation 



*fc:0 



ip-i(k)al :1 + t/>-2(fc)°*:2 + i>-a(k)al. 3 , 

ipoi(k) a\. x + ipo2(k) ak-.2 + ^03 (fc) ajU, ( 4 ) 



a[. + = ip + i(k)a k :i +ip +2 {k)a\. 2 + ip +3 (k)a k , 3 , 
we reach a spin-wave Hamiltonian, 

H = E g + ^2 v\(k)al;\Oik:\, 



(5) 
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FIG. 2: Dispersion relations of the elementary excitations 
in the spin-i plaquette chain, two (circles and diamonds) of 
which reduce the ground-state magnetization and are thus 
of ferromagnetic character, while the rest (squares) of which 
enhances the ground-state magnetization and is thus of anti- 
ferromagnetic character. The exact-diagonalization results at 
N = 4, 6, and 8 are presented by symbols of small, middle, 
and large sizes, respectively, whereas the up-to-0(S' 1 ) linear 
(LSW) and up-to-0(5' ) interacting (ISW) spin-wave calcu- 
lations are given by dotted and solid lines, respectively. 
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FIG. 3: (Color online) Dispersion relations of the elementary excitations in the spin-i trimeric chain with varying 5. The 
first-order perturbational calculations, together with the up-to-0(5 rl ) linear spin-wave findings, are given in the upper five, 
whereas the second-order perturbational calculations, together with the up-to-0(S' ) interacting spin-wave findings, are given 
in the lower five. The exact-diagonalization results at N — 4, 6, and 8 are presented in both upper and lower panels by symbols 
of small, middle, and large sizes, respectively 
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With E g --: J2i=2,l,0 E S and ^A(fe) = Ei=l,0 W A ( fe )> 
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state energy, while Eg' and u> A (fc) (i = 1, 0, • • •) are the 
0(S l ) quantum corrections to the ground-state energy 
and the dispersion relation of mode A, respectively. Here 
we have discarded the 0(S^ 1 ) terms. There are sev- 
eral way s 40 i 41 ' 42 i 43 ' 44 of treating the quartic interactions. 
When we diagonalize the one-body terms and then take 
account of the two-body terms perturbationally,— the 
spin-wave energies read 
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and their eigenvectors are given by 
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Figure [2] shows the thus-calculated spin-wave ex- 
citation modes together with the exact eigenvalues. 
Free spin waves well describe the ferromagnetic modes 



uj-(k) and u>o(k), while higher-order quantum correc- 
tions play an essential role in reproducing the antifer- 
romagnetic mode u/+(k). The 0(S ) quantum correc- 
tions significantly improve fully delocalized magnetic ex- 
citations in general ) 18 ' 41 i 46 but the standard Holstein- 
Primakoff magnon series expansion seems not to work 
well for highly localized excitations. The dispersive 
branches w^(fc) are nothing but the elementary excita- 
tion modes of spin-alternating linear-chain Heisenberg 
ferrimagnets.— They are parallel within the spin-wave 
theory, but their difference w+(fc) — w_(fc) is momentum 
dependent in fact. On the other hand, the flat band 
uio(k) arises from further excitation degrees of freedom 
in the present system. When J\ — J2, the Hamiltonian 
([1]) reads 



N 



posite spins Tn:3;n+i:i- All the composite spins are sat- 
urated in the ground state, r n:3 . n+1:1 = 25(25 + 1), 
and therefore, their excitations are necessarily of ferro- 
magnetic aspect. The ferromagnetic excitations of local 
character are well understandable within the spin-wave 
description. Equations (@|) and (TlO|) show that a' n:1 and 
creating bosonic excitations on the Cu 2+ (2) sites, 



l n:3i 



indeed participate in the construction of o:I. , but any of 






n:3;n+l:l , 



(13) 



creating bosonic excitations on the Cu 2+ (1) sites, 
does not. Without mediation of bridging spins S n: 2, any 
intraplaquette excitation is never movable. Then what 
may happen with S moving away from unity? The spin- 
wave theory, whether up to 0(S 1 ) or up to 0(S°), pre- 
dicts that the ferromagnetic and antiferromagnctic exci- 
tation modes w^(fc) are still parallel and the extra ferro- 
magnetic excitation mode between them, wo(fc), remains 
dispersionless. Let us verify the true scenario. 



with composite spins r w :3;ri+i:i — 5„:3 + Sn+i-.i, each 
lying diagonally across an elementary palquette. Since 
the Hamiltonian ([TBI) commutes with T 2 , 



n:3;n+l:l 

T n :3;n+i:i(Tn:3;n+i:i + 1), we have good quantum num- 
bers r n :3; n +i : i, each taking cither or 1. Therefore, the 
plaquette-chain Hamiltonian is block-diagonalized by the 
set of numbers {T n: 3 ;n _|_i : i; n = 1, 2, • • • , TV} 4 ^ as well as 

by the total magnetization J2n=i( S n-.2 + T n:3 ;n +i:i) = 
M. The Hilbert space of En=i( T n : 3;n+ia) 2 /2 = 
E„ = i(*n:3 • 5 n+ i:i + 3/4) = JV = N corresponds to 
the ferrimagnetic chain of alternating spins 1 and \ 
and consequently we here have exactly the same disper- 
sion relations 4 ^ of elementary excitations. The Hilbert 
space of M — N — 1 and M. — N/2 — 1 consists 
of N subspaces labeled {Ti :3;2: i, 12:3:3:1, •••■> Tn:3-,i-.i} = 
{0,1,---,1},{1,0,1,---,1}, ■■■,{!,■■■, 1,0}, and they 
all give the same set of eigenvalues, forming N flat bands. 
We find the lowest one in Fig. O 

Thus and Thus, the spin-5 plaquette chain turns out a 
combination of the alternating-spin- (25, 5) linear chain 
and extra excitation degrees of freedom within the com- 
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FIG. 4: Probability of two spin i's constructing a spin 1 in 



the ground state of the spin-i trimeric chain of N — 64 with 
varying 6, where P„:3;n+i:i = T^ :3;n+1:1 /2 = S n -.3 ■ S n+ ui + 
3/4 and P n: i;n;3 = T^ :1 . n:3 /2 = S n -.i -5 n:3 + 3/4 are estimated 
by a quantum Monte Carlo method. 



III. BOND-ALTERNATING TRIMERIC 
CHAINS 

We demonstrate in Fig. [3] several schemes of calcu- 
lating low- lying excitation modes for the spin-i bond- 
alternating trimeric chain. In spite of the persistent flat 
band within the spin-wave theory, the exact diagonal- 
ization reveals that it can be dispersive with varying S. 
When (5^1, the Hamiltonian ([I]) does not commute 
with T n . 3 . n+l . l . Now that there is a certain probability 
of composite spins T n: 3 ;n +i:i being singlet even in the 
ground state, the gapped ferromagnetic excitation mode 
u>o(k) is no more describable as their individual triplet- 
to-singlet flips. At 6 = 0, any excitation is localized 
within a trimcr of S n: i, S n -2, and 5„ : 3, and the excitation 
spectrum degenerates into three flat bands, w_(fc) = 0, 
u>o(k) — Ji, and 0J+(k) = 3Ji/2. Figure [3] shows that 
the middle branch of them connects with the flat band 
at S = 1. Without J2, the Hamiltonian (fl]) is reduced to 



JY 



T~t — Ti-i — J\ 2_^ 5 ra: 2 



n:l;n:3, 



(14) 



71=1 



with intratrimer composite spins T„ : i ; „ : 3 = S n -.i + 5„ : 3 
and thus commutes with T n:1 . n:3 . The plaquette chain 
(fT3"]) and the decoupled trimers l[14"|) both exhibit a flat 
band due to gapped ferromagnetic excitations, but their 
ways of constructing local immobile excitations are dif- 
ferent from each other. Triplet-to-singlet [(45 + l)-fold- 
multiplet-breaking in general] flips of intraplaquette com- 
posite spins T n: 3 ;n _|_i : i are the elementary excitations in 
the former, while those of intratrimer composite spins 
T n: i;n:3 are the elementary excitations in the latter. Fig- 
ure [3] shows how such composite spins behave in the 
ground state with varying 6. Neither T„ : 3 ;n+ i : i nor 
T n :\- n :3 form complete triplets at < 6 < 1, due to 
nonvanishing off-diagonal matrix elements (T„ : 3 ; „_|_i : i = 
l|W|T n:3i „+i:i = 0) and (T„ :1; „ :3 = l|W|T n:1;n:3 = 0). 



At the two particular points 5 — 1 and 6 = 0, only 
the Cu 2+ (2) ion spins S n: i and S n:3 constitute the 
gapped ferromagnetic excitation mode, but otherwise the 
Cu 2+ (1) ion spins S n -2 also contribute to that. Without 
interconnecting spins S n: 2, any excitation is immobile, 
whereas with their mediation, all the local excitations 
can be itinerant and the resultant bands are dispersive. 
Perturbational calculations support such a scenario. 
With increasing couplings J 2 between isolated trimers, 
the energy dispersion relations grow as follows: 



E S 1 6 869 ,2 ^/riJN 

^ = ^(1 -cos fc) + 4(929-474 



x cos k - 455 cos 2k) + 0(S 3 



(15) 



(16) 



Ji 



1 - 0.38970(5 + <5 2 (0.32099 - 0.26736 



x cos k + 0.04212 cos2fc) + 0(6 3 ), 



(17) 



^) = 3 + ^ (7 _ 12co8fe) + ^ (346 _ 10 

xcosfc- 109cos2fc) + O(<5 3 ). (18) 

Equations (|16[) - (f]~8|) are also drawn in Fig. [3] The first- 
order perturbation points out that not only ui^(k) them- 
selves but also their difference should be dispersive, but 
it cannot reveal nonvanishing momentum dependence of 
u>o(k). We cannot reproduce the dispersive middle band 
until we take account of the second-order perturbation. 
The lowest ferromagnetic and antifcrromagnetic excita- 
tions of decoupled trimers (fT4"| . gapless and gapped by 
3Ji/2 from the ground state, respectively, are both N- 
fold degenerate and are expressed as 

|£; t (to)) = |-(1±3)J 1 /4;1/2t1> 



n^m 



|-Ji;l/2) B (m=l,2,..-JV), (19) 



while their gapped ferromagnetic excitations at an energy 
cost of J\ are 7V 2 -fold degenerate and are expressed as 

\E (m,m')) = 8 mm ,\0;-l/2) m <g> |-J i; l/2)„ 

+ (1 - S mm ,)\0- l/2) m O | - Ji; -l/2) m » 

® |-Ji;l/2)„ (m,m' = 1,2, •••JV), (20) 

in terms of the eigenstates of an isolated trimer 
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(|TU) + |m> + IUT)), 
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With perturbational interactions 7^2 turned on, the 
iV-fold degeneracy of the eigenvalue —[N — 3(1 =F 
l)/4]Ji = (E T (m)\Hi\E T (m)) is completely lifted, 
whereas the iV 2 -fold degenerate eigenvalue —(N— 1) J\ = 
(E (m,m')\Hi\Eo(m,m')) only splits into N flat bands 
within the first-order corrections. The second-order cor- 
rections are necessary for reproducing the dispersion re- 
lation of u>o(k). In this context we may be reminded that 
Honecker and Lauchli 49 pionceringly investigated anal- 
ogous but frustrated Cu 2+ trimeric chains. The gap- 
less ferromagnetic excitation mode (J16I) is indeed derived 
from their effective Hamiltonian under strong trimeriza- 
tion 6 < 1. 



IV. SUMMARY AND DISCUSSION 

We have investigated the low-energy structure of inter- 
twining double-chain ferrimagnets composed of homoge- 
neous spins with particular emphasis on the gapped fer- 
romagnetic excitation mode. While there exist a macro- 
scopic number of flat bands^ in the excitation spectrum 
at 6 = 1 and 6 = 0, which signify uncorrelated excita- 
tions of local spin-2S' multiplets in any case, they become 
dispersive as soon as 6 moves away from these particular 
points. Pair excitations of corner spins S n: 3 and S n +i:i 
are elementary in plaquette chains of 6 = 1 , while those 
of S n: i and S n .3 are elementary in decoupled trimers of 
(5 = 0, both of which are completely immobile without 
any mediation of joint spins S n -.2- The spin- wave theory 
successfully characterizes the plaquette chain but fails 
to find arising contribution of S n: 2 to gapped ferromag- 
netic excitations with bond alternation. Such a mislead- 
ing prediction has been corrected by further numerical 
and analytical investigations. 

The spin-5 plaquette chain thus shares whole the na- 
ture of the altcrnating-spin-(25, S) linear chain and fur- 
ther exhibits ferromagnetic excitations of its own. All 
the findings but the flat band in Fig. [5] are indeed ex- 
actly the same as we have in the ferrimagnetic Heisenberg 
chain of alternating spins 1 and ^. 45 Even though the 
homogeneous-spin plaquette chain and the alternating- 
spin linear chain are equivalent in their ground states, the 
former demonstrates its extra excitation degrees of free- 
dom and deviates from the latter with increasing temper- 
ature. In order to illuminate similarities and differences 
between them, we show in Fig. [5] quantum Monte Carlo 
calculations of their static structure factors 



n,l,n',l' 



ig(x n ..i—3 



) qz qz 



(22) 




W 



0.0 
0.0 



A a 7V, 



m 


■ 


i i ■ 


so 




O 


.5.0 


°x x x° 



0." 



1,0 



II 



k a TtJ, - 0. 1 a k B TlJj - 0.3 n k B TlJi - 0.6 
x k a T/Ji ■ 0.2 k B TU\ ■ 0.4 - k B TUi "1.0 





10 2.0 

FIG. 5: (Color online) Quantum Monte Carlo calculations of the static structure factor S(q), with the distance between 
neighboring spins in the chain direction set equal to unity, as a function of temperature. The whole view and enlargements at 
q = and q = n for the spin-i plaquette chain of TV = 64 (a) and the alternating-spin-(l, |) linear chain oi N = 64 (b). The 
ferromagnetic [5(0)] and antiferromagnetic [S(7r)] peaks are observed in more detail (c), where the common asymptotic values 
in the high-temperature limit, 3/4 and 11/12 for the the spin-i plaquette chain and the alternating-spin-(l, ~) linear chain, 
respectively, are indicated with arrows. Thermal averages of the projection P n -.3;n+i:i = T^. 3 . n+1 . 1 /2 = S n -.3 ■ S n +i-.i + 3/4 
in the spin-^ plaquette chain are also shown for reference, where the asymptotic value in the high-temperature limit, 3/4, is 
indicated with an arrow. 



as functions of temperature, where the chain-directional 
coordinates x n .j are given in the unit of neighboring-spin 
spacing. The pronounced peaks at q — and q = tt 
reflect the ferromagnetic and antiferromagnetic double 
excitation mechanism in common. Without any field 
applied, 5(0) and 5(7r) are, respectively, the uniform 
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FIG. 6: (Color online) Quantum Monte Carlo calculations of 
the static structure factor S(q), with the distance between 
neighboring spins in the chain direction set equal to unity, as 
a function of bond alternation and temperature for the spin-i 
trimeric chain of N = 64. 



and the staggered susceptibilities multiplied by temper- 
ature. With decreasing temperature, they both diverge 
as 1/Ti 38 i 45 i 51 With increasing temperature, they both 
approach the paramagnetic value J2i S n :i(S n -.i + l)/3 but 
behave differently at intermediate temperatures. A mini- 
mum of 5(0) as a function of temperature is characteristic 
of ferrimagnetsi£ii&2iii2{ii££i^i£i 5(0) monotonically de- 
creases and increases with increasing temperature in fer- 
romagnets and antiferromagnets, respectively^ Though 
the thermal as well as quantum behaviors of the spin- 
5 plaquette chain and the alternating-spin- (25, 5) lin- 
ear chain are very much alike, yet there grows a differ- 
ence between them with pair excitations of intraplaque- 
tte spins 5„ : 3 and S n +i:i from their highest multiplets. 
The ferromagnetic and antiferromagnetic structures of 
S(q) less survive increasing temperature in the spin-5 
plaquette chain than in the alternating-spin-(25, 5) lin- 
ear chain. The larger 5 the major difference in 5(g) as 
limT-oo^^fe) - S (s ' s < s \q)} = 25 2 /3. Alternating- 
spin- (25, 5) ferrimagnetic chains behave like combina- 
tions of spin-5 ferromagnetic and spin- (25) antiferromag- 
netic chains^ while such a simple magnetic sum rule 
is not available to intertwining double-chain ferrimag- 
nets of our interest. Additional intraplaquette antifer- 
romagnetic interactions induce incommensurate peaks in 
5(g), 55 making corner spins S n: s and S n +ia frustrated. 

Once S moves away from unity, the homogeneous- 
spin trimeric chain never more shares any feature of the 



alternating-spin chain. Figure [6] presents S(q) with vary- 
ing 8 and analyzes its features at q — and q = tt in 
particular. At low temperatures, 5*(0) and S(ir) both 
decline with decreasing 8, but they still diverge as \/T 
unless 8 = 0i22 At high temperatures, S(ir) remains de- 
creasing, whereas 5(0) turns increasing, with decreasing 
8. Decoupled trimers are nothing more than paramag- 
nets and their structure factor is given as 



>Jl/k B T 



-J 1 /2k B T 



qM = 3 _ 2 

{q> 4 3e J i/ fc BT + i + 2e-- 7 i/ 2fe BT 

1 e Ji/k B T _ 3 + 2e -J 1 /2k B T 



cosq 



strong bond alternation 8 < q. 1 ,30,32 ,33^4,35 it may be 
hard to detect the dispersion relation ujo(k) there. Ni 
analogs, if available, will present an energy structure of 
the same type on an enlarged energy scale. Highly lo- 
calized excitations in fully exchange-coupled bulk mag- 
nets may either arise from an accidental arrangement of 
exchange couplings or come out of a particular lattice 
structure of geometric aspect. While some examples of 
the former case can be found in spin-i bond-polymerized 
chains in principle ; 3 ^ 61 it must be hard to observe them 
in real materials. On the other hand, the present find- 
ings are of the latter origin and may be more accessi- 



6 P Ji/k B T I i I o 9 -Ji/ik B T cos ^' (^) ble experimentally. The Shastry-Sutherland lattice 5 ^ is 



which is also drawn in Fig. [S] with solid lines. Equa- 
tion (|23| at q — reads as the effective Curie law for a 
trimer entity, where the Curie constant varies from 1/4, 
attributable to the ground-state doublet, to 3/4, simply 
coming from free spin i's, with increasing temperature. 
Arising intertrimer couplings J2 immediately pronounce 
a quadratic dispersion relation of the ferromagnetic exci- 
tations at small momenta and their further increase costs 
the antiferromagnetic excitations higher energy. That is 
why growing global correlations enhance and reduce the 
uniform susceptibility-temperature product at low and 
high temperatures, respectively. 

Weak but nonvanishing dispersion of the gapped fer- 
romagnetic excitation mode is the most remarkable find- 
ings of ours and is the very characteristic of intertwining 
double-chain ferrimagnets. As the existent compounds 
^3Cu3(P04)4 have all been reported to exhibit rather 



also interesting in this sense and the model compound 
SrCu2 (1303)2 indeed exhibits an excitation mode of lit- 
tle dispersion in its low-energy spectrum. 57 Such local 
excitations crystallize to form a superlattice, quantizing 
the ground-state magnetization i 58 i 59 ' 60 There may be a 
similar scenario in low-dimensional ferrimagnets of topo- 
logical origin as well. We hope the present study will 
stimulate further experimental explorations. 
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